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We have systematically studied the S-wave pseudoscalar meson and heavy vector meson scattering 
lengths to the third order with the chiral perturbation theory, which will be helpful to reveal their 
strong interaction. For comparison, we have presented the numerical results of the scattering lengths 

■ (1) in the framework of the heavy meson chiral perturbation theory and (2) in the framework of the 
' infrared regularization. The chiral expansion converges well in some channels. 
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^ I I. INTRODUCTION 

In the past eight years we have witnessed the renaissance of the hadron spectroscopy. Many interesting new hadron 
states were discovered experimentally, some of which do not fit into the quark model easily. These new hadron states 
include (1) the famous XYZ states, which are either charmonium or charmonium-like states above the open-charm 
decay threshold; (2) the narrow charm-strange mesons D s o(2317), D s i(2460) etc.; (3) the charged Upsilon-like Z^ 
\ states recently announced by Belle collaboration A general feature of many of these new hadrons is that they 
are very close to the two open- flavor meson threshold. For example, X(3872) is very close to the DD*, pJ/ip and 
ujJ/ip threshold. _D s o(2317) is very close to the DK threshold. These charged Z^ states are very close to the B^B* 
threshold. Because of their proximity to the two meson threshold, one may wonder whether some of these new hadron 
states are good candidates of loosely bound molecular states composed of two mesons? Or will the coupled-channel 
effect between the bare qq state in the quark model and the two meson continuum help lower and push the mass of 
the bare quark model state close to the threshold? In order to allow the above two mechanisms work, there must 
exist attractive interaction between the two mesons. 

Generally speaking, it is very difficult to study the hadron interaction starting from the first principle of strong 
interaction, i.e., quantum chromodynamics. Most of such investigations are performed on the lattice numerically. 
However, we may turn to chiral perturbation theory for help if one of the interacting mesons is a light pseudoscalar 
meson. In this case one can derive the scattering amplitude order by order rigourously. From the scattering amplitude, 
we can extract the scattering length, which is directly related to the hadron interaction. In this work we will study the 
pseudoscalar meson <fi and heavy vector meson D* scattering lengths in order to learn whether there exits attraction 
between the pseudoscalar meson (f> and heavy vector meson D*. Such a study will provide valuable information on 
their interaction to the D*K system. 

Up to now, a few new charmed mesons and their properties have attracted much interest over the past few years, 
especially the extremely narrow Z? s o(2317), Z? s i(2460) states. The D(D S ) and D*(D*) mesons constitute the lightest 
charmed doublet according to heavy quark symmetry. The recently observed new charmed particle -D s o(2317) is 
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speculated to be a candidate of possible molecular states composed of the D meson and K meson. The DK interaction 
is very important for us to understand the underlying structure of the D s o(2317) meson 0, There have been some 
lattice investigations on the ttD scattering 

[4 6], KD scattering and D s o(2317) |7H9|. If there are strongly attractive 
interactions between them, the D* meson and pseudoscalar meson might also form possible molecular states. Our 
present study of the <pD* scattering with chiral perturbation theory would be helpful to the future lattice simulation 
of the <j>D* scattering numerically. 

Chiral perturbation theory and lattice QCD are widely used to study the hadron interaction in the nonperturbative 
region of QCD [IM3- 

The interaction between the D meson and the light pseudoscalar meson was studied recently 
with chiral perturbation theory p^| - |20j . It is interesting to extend the same formalism to study the interaction of 
the heavy vector meson D* and the light pseudoscalar meson. In this paper, we will calculate the S-wave scattering 
lengths of the pseudoscalar meson and D* meson with the heavy meson chiral perturbation theory (HMxPT) and 
the infrared regularization f IR) [2l| - [28j . The scattering length a^D* reflects their interaction. In our formalism a^w 
is related to the threshold T-matrix T^d* ■ T^d* — 8tt(1 + )a<j>D* , where and Mo* are the masses of the light 
pseudoscalar meson and D* meson respectively. 

With the explicit power counting scheme, rlM^PT is a useful tool to investigate the heavy meson interactions 
plM^ . We will expand our calculation by e = p/A x , where p represents the momentum of the light pseudoscalar 
meson, the small residue momentum of the heavy meson in the nonrelativistic limit, or the mass difference 5 between 
D and D* mesons, and A x represents either the chiral symmetry breaking scale around Awf^ or the heavy mesons' 
masses M (about 1900 MeV) in the chiral and heavy quark symmetry limit. The IR scheme is also a useful tool 
based on chiral perturbation theory [25l - |28j , which ensures both good power counting and correct analyticity. IR and 
HMxPT generally lead to the same results except that the IR formalism includes the higher-order infrared parts of 
the loop graphs (251 ] . 

This paper is organized as follows. In Sec. II we list the basic notations, the relevant Lagrangians, and the chiral 
corrections to the threshold T-matrices with HMxPT. We present the IR expressions in Sec. III. The low-energy 
constants (LECs) are estimated in Sec. IV. Finally, we give the numerical results and discussions in Sec. V. 



II. THE T-MATRICES WITH THE HEAVY MESON CHIRAL PERTURBATION THEORY 



We list the Lagrangian of HMxPT at the leading order here, 



& = / 2 Tr («/ + *±) , (1) 



,(i) .,. .„ ^ .„ , ^ 1 



£n* = ~(( iv ■ dH )H) + (Hv ■ TH) + g(H faB) - ^S(Ha^Ha^), (2) 
where / is the decay constant of the pseudoscalar meson in the chiral limit, and 

H = {P*^ + iPjs) , H = 7 °iiV = (p; V + iPS) ^> (3) 

P=(Lfi,D+,D+), P; = (D*°,U*+ (4) 

Heavy quark symmetry is exact only when the heavy quark mass is infinite. In this work we will also systematically 
include effects of the explicitly broken heavy quark symmetry through the last term containing the D* and D mass 
difference 5 in Eq. @. The notations read 

1^ = ^,^], u,= l -{^,8^}, £ = exp(#/2/), x±=Z j X?±txt, X = ^g(rnl,rnl,2m%-ml), (5) 



= V2 n -^ + ^ K° . (6) 




The following Lagrangians at the second and third order are used in the calculation of the threshold T-matrices, 1 

,(2) 

-H<f> 



£<£l = c (HH)Tr(x + ) + c 1 (Hx + H)-c 2 (HH)Tr(vuvu)-c 3 (HvuvuH), (7) 



1 The sign in front of ci in Eq. (9) of Ref. 18] should be +. The signs in Eqs. (12) and (13) should be consequently changed. 
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C^l = K 8{HH)Tr(x+) + K 1 S(Hx+H)-K2S(HH}Tr(vuvu)-K 3 S{HvuvuH) + K(H[x-,v-u]H). (8) 

The 0(e 2 ) and 0(e 3 ) Lagrangians could also contain terms like (Ha^ H<t^)Tt(x+), (Ha^x+Ha^). These terms 
break the heavy quark symmetry hence are suppressed. They lead to different LECs c^'s and k^s for the <j)D and <j)D* 
scattering, although they do not result in the new independent vertices we need. 

There are eleven independent T-matrices in the pseudoscalar meson and D* meson scattering due to the isospin 
symmetry. The threshold T-matrices start at 0(e), which can be derived from Eq. ([2]) 



T (3/2) 
± ttD* ~ 




T (l/2) 
1 ttD* 


2m v 


- t 7r£>; 


= o, 


T (0) 

1 KD* — 


2niK 


T (l) n 
1 KD* — u > 


T (l/2) 
± KD* — 


m K 


(2 ' 

TT 


(2 ' 

■1 7T 


f2 ' 
J K 


' Pk 


T (l) 
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m K 
Pk' 


T (o) 

KD* 


m K 

- f2 ' 
J K 


r (l/2) 
KD* 


m K 


T (l/2 
' 1 tjD* 


= o, 


T (0) _ Q 







(9) 

where the superscript in the bracket represents the total isospin of the channel. We express T-matrices with the 
renormalized decay constants f^, fx and f v [1^, [3(| rather than / here. The difference could be accounted for 
through T-matrices at 0(e 3 ) or higher order. 
Similarly we get the results at 0(e 2 ), 2 

T (3/2) _ 8c to 2 + 4cito 2 + 2c 2 to 2 + c 3 to 2 T (i/2) _ 8c to 2 + 4cito 2 + 2c 2 m 2 + c 3 m 2 

ttD* f2 ' ttD* r2 ' 

J TT J TT 

(!) _ 8c m 2 + 2c 2 m 2 (0) _ 8c m 2 K + 8cim 2 K + 2c 2 m 2 K + 2c 3 m 2 K (1) _ 8c m|- + 2c 2 m|- 

ttD* ~ f2 ' KD* ~ f2 ' KD* ~ f2 ' 

Jtt Jk JK 

(1/2) _ 8c m 2 K + 4ci m|- + 2c 2 m|- + c 3 m|- ( X) _ 8c m|- + 4gm 2 K + 2c 2 m 2 K + c 3 m 2 K 

1 KD* S ~ f2 ' KD* ~ f2 ' 

•/if JK 

(0) _ 8c m|- - Aciui 2 ^ + 2c 2 m|- - c^m 2 ^ T (i/2) _ 8c m 2 K + 401771^ + 2c 2 m|- + c 3 w|- 

— (2 ' J i?D* — f2 ' 



,(1/2) _ 24c m 2 + 4cim 2 + 6c 2 m 2 + c 3 m 2 (0) _ 24c m 2 + 32cim|- - 16cim 2 + 6c 2 m 2 + 4c 3 m 2 

rjD* ~ 3 J2 ' ^-D; — 3 y2 ' 



Here we have used the Gell-Mann-Okubo mass relation m 2 = {Am 2 K — m 2 )/3, which makes the expressions more 
concise. 

The T-matrices contain contributions from both the tree and loop diagrams. We show all eighteen loop diagrams 
which contribute to the threshold T-matrix at 0(e 3 ) in Fig. [TJ We calculate them with the dimensional regularization 
and modified minimal subtraction. More specifically, for the unstable D* meson we renormalize its wave function at 
the point f 



d(r) 



(11) 



where Ii£)*(r) is the one-particle irreducible D* self-energy, r is the remaining momentum r = v ■ p — M, and f is 
the complex pole of the propagator f — 5 — Hd* (r) = 0. The divergences from loops can be absorbed after the wave 
function renormalization and redefinitions of Ki, 



2g 2 L A . r 5g 2 L „ 3g 2 L 3T 

4:Ko + K 2 = -g^2-+4K + K 2 , K 1 = Y^2 +K 1' K3 = 'jT +K 3> K= TJ2 +K ' ( l2 > 



where 



xD-4 



1 



L = -{ + -h E -l-ln4ir)} . (Euler constant lE = 0.5772157) (13) 

l07r 2 I D — 4 2 J 

Here A is the scale of the dimensional regularization. We will set it at 4^/^, AttJk and 471"/,, respectively for the pion-, 
kaon- and 77-scattering. 



2 The constant Ci in of Ref. [Ill should be ^(3Ci — Co). The corrected T-matrix is the same as T?|j» here. 
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In order to make the expressions short, we introduce the following notations and functions: 



J = ~ 247r2 {m J_ mJ) /4 ( W + W - 6^ cos-H-^) + 3^ log ^ - 2^ 5 

_ 2< 5 3 log ^ - 2nm n 3 - 2(5 2 - m^f' 2 log ^ + 2(<S 2 - m 2 fl 2 log V ^ ~ ^ + - 
A A A 

TO,, 



-2i7r5V<5 2 - Wtt 2 + 2i 7 rm 7r 2 V'5 2 -m 7r 2 - Sm^J log + 2m 7r 2 5 + 2<5 3 log 

A A 

g 2 (2Vm 2 -S 2 cos- 1 (-£) +2(5 log (f) + 2ttto - <5 m > 5 

= ~16tt 2 / 4 \2<Jlog f + 2V^ 2 - m 2 (log VA "-^ +<5 - iir) + 2vrm - J m < S ' 

3i m 3 2 f- Vra 2 -U > 2 COS- X (-^) TO 2 > OJ 2 

V(m,u) = - -4- - 4- j V^^log m 2 < w 2 , W <0. (14) 



tt 2 / 4 2tt 2 / 4 n 2 f 



VCJ 2 - m 2 (- log V^-m^+ H + ^ to 2 <w 2 jW >q 



' 1 



FIG. 1: Nonvanishing loop diagrams for the pseudoscalar meson and D* meson scattering lengths to 0(e 3 ) with HM^PT and 
IR. The dashed lines, thin solid lines and thick solid lines represent the pseudoscalar Goldstone bosons, D mesons and D* 
mesons, respectively. 



In the third-order T-matrices there are some terms proportional to those in Eq. (|10[) . Therefore, we divide the 
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T-matrices into two parts, 

T = T + T 2 , (15) 
where T 2 can be obtained from Eq. (jlOp through the following replacement: 

c o K o$, c i k\5, C2 — > k 2 5, C3 — > k^5. (16) 

The remaining T at 0(e 3 ) reads 

Tin? = -\v{m K , -m n ) - lv(m n , -m w ) - \v(m*, m,) + ^-mlW{m v ) - £m= W(m*) + 2V 1 - ^§^", 

8 8 8 18 I 

1 3 111 16 ^ ^ 

= -^V{m K -m n ) - —V(m K ,m n ) - -V(m n ,m w ) + -m^(m^) - -m^(m w ) - AV l + ™f , 

1 12 

T^, = --V(mjf, -m,) - gV(mj<-, m,) + -mlw(m n ), 



3 11 1 16m 3 K r 

-V(m n ,m K ) - -V(m K ,-m K ) - -V(m K ,m K ) - Jm 2 K - ^m K W{m v ) - W 2 - 3V 3 + V 4 H 7^ 

8 8^ y j K 



f2 



2^],. = -in^,™K)-^(m K ,-™ K ) + ^-i77i 2 K ^(m r) )-^ + y 3 + y, 



./K 

3 3 11 1 

Trd* = 7^ V ( m n> - m K) - TT-yim^, -m K ) + -V(m K , -m K ) - -V{m K ,m K ) + Jm 2 K - -m 2 K W{m v ) - 3V S 



16 v JW 16 v "' 4 v ' 8 v " y ^9 

J Is 
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~*Rd* = ~J^ v ( m rn m K) - j^V(m n ,m K ) - ^V(m K ,-m K ) - ^V(m K ,m K ) - ^m 2 K W(m tl ) - W 2 + 8m ^ K , 



^d* 2) = -^ V '( m ^' _m ')) ~ Ye^™^'™ 1 ?) _ ^ m2 K W ( m K) + -m 2 W{m v ) - —mlW(m v ) + -mlW(m„), 



f i% = -^n^,-™r,)-^n^ ! ^)-^K^(^) + ^m^(m,)-^m2iyK), (17) 

where Vj. = V 2 = V 3 = V4 = with HM^PT but nonzero with IR. 

The T-matrices of the pseudoscalar meson <\> and D meson scattering with HMxPT were derived in Ref. [18| ■ For 
the 4>D case, there is no </>_DD-vertex in the leading order. Therefore in the loop calculation there are no similar 
diagrams of the third column in Fig. [1] where <f>D are the intermediate states. If we let 5 — > and neglect the explicit 
breaking of heavy quark symmetry, the T-matrices of the pseudoscalar meson and D* scattering will be the same 
as those of the pseudoscalar meson and D meson scattering to 0(e 3 ) at the threshold [l|j], which is required by the 
heavy quark symmetry. 

III. THE T-MATRICES WITH THE INFRARED REGULARIZATION METHOD 

For the IR scheme, we use the heavy meson Lagrangian with the relativistic form at the leading order, 

L h\ = V^PV^P^ - M 2 PP^ - V^P^V^Py + (M + 5fP*"P? 

+z2. 9 M(P;^pt _ p u »p^) + g(P*u a VpP? - VpP;u a P?)e» vafi , (18) 
where heavy quark symmetry is also assumed to relate the couplings of the nD* D- vertex and 7rD*D*-vertex, and 
1 ;-n i> — :;» i> 1 '""/>. ;n i>\ — ;;i i>\ 1 1 



P = -f^, W tl P a =id ll P a -T b «P by mpP*=idpf*+r?Pl (Similar for P*.) (19) 



M 
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The second-order and third-order Lagrangians we use are the same as the relevant terms in Eqs. ([Jj [8]), but some 
coefficients should be redefined to fit the experimental data, 

C%\ = 2M (c P*' 1 P; t Tr(x + ) + dP*"^ - c 2 P^P^Tr(u ■ u) - c 3 P*»u ■ uP?), (20) 
4g = 2M(KoP^P^Tr( X+ )6 + KtP^x+P^S - K 2 P^P^Tv(u ■ u)S - k 3 P*»u ■ uP?S) 

+Mk {v v P*»[x-,u v ]Py - P*»[x-,u»]V v Py) . (21) 

The T-matrices are nearly the same as those of HMxPT except that the expressions of J, W(m), V(m, uj), V\, V2, 
V3, and V4 are more complicated. We list their definitions in the infrared scheme in Eq. (|B7[) in the Appendix. We 
have also verified that the results with IR are the same as those with HMxPT when M approaches to infinity. 



IV. LOW-ENERGY CONSTANTS 



The difference of our results between HMxPT and IR originates at the third order due to different loop integrals. 
We use the same LECs at the first and second order for both HMxPT and IR. At the leading order, we have |31| 

ni-x = 139 MeV, m K = 494 MeV, 5 = 142 MeV, 

/ 7r = 92MeV, Jk = 113 MeV, f v = 1.2f K , 9 = 0.59. 

At the second order, from the mass splitting between heavy mesons we get 3 

ci = — — ^ ^ -- = 0.12GcV- 1 . (22) 

16 M(m 2 K -ml) 

In order to obtain other LECs that can not be determined from the available experimental data, we resort to the 
resonance saturation model (32l . l33j . 

At 0(e 2 ) only the light unflavored mesons with J p = + and charmed mesons with J p = 1 + con- 
tribute to the 4>(f)D*D* vertex at threshold. Thus we consider the scalar singlet er(cr(600)), the scalar octet 
k(k(800), ao(980), /o(980)), and Z? s i(2460) vector triplet in this section. In the Appendix [A] we will discuss the 
uncertainty of the LECs at this order. 

Here we list the corresponding effective Lagrangians, 

Ca-n-K = 4c rf Tr(u • m)ct + c m Tr(x+)fx, 
C aP »p« = c a P*»P^a. (23) 

Integrating the a meson out through the t-channel we get 

^ ^Tr( x+)W; f + ^Tr(v -uv u)P"*P? (24) 

Similarly, from the Lagrangians of the scalar octet n 

Cktvtt = 4c d Tr(u • uk) + c m Tr(x+K), 
4p.p. = c k P*^kP;\ (25) 

one obtains 

C K eS ~ -^Tr( X+ )P*»P? + ^P*»X + P? ~ ?p?Tr(t, -uv u)P*»P? + ^P^v -uv uP? (26) 
Integrating D s i(2460) out of the following interacting Lagrangian, 

£o sl (2460) = Gi (D^(2m)u v id u P; - id u P^u v D^(24Q0)) + G 2 (-id v D^{2A60)u u P; + P;« w ift / £)£i(2460)) 



3 The coefficient 4 in Eq. (33) of Ref. [3 should be 8. The correct values of ci are 0.12 GeV -1 in Eq. J22J and 0.10 GeV -1 in Eq. J36) 
of this paper. 
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= (G 1 +G 2 )(D^(2460)u^P;~ l ^P; t u^ 1 (2460))+O(e 2 ), (27) 

one gets 

= - \? 2 X+G2?M °l P»>v -UV uPf. (28) 

IV1 D sl (2460) ~~ 1X1 D* 

The effective coupiing constants \G\ + G 2 | were estimated with QCD sum rule approach in Ref. [34]: \G\ + G2I = 
1.2 ±0.2. 

Adding the above effective Lagrangians (121)1) and (|2"5)l together, one can estimate the LECs by comparing the 
sum with the relevant terms in Eq. ([7]), 

Co- Cm c K c m c K c m 

C2 = 4Cj5d + 4c " Crf C3 = 4c * Crf + |Gi + G 2 | 2 M D » 
For the broad resonances er(600) and k(800), we use the masses and widths extracted from a model-independent 

way m m, 

771 (j = 44llg 6 MeV, r CT = 544+^ MeV; m K(800) = 658 ± 13 MeV, r K(800) = 557 ± 24 MeV. (30) 

In our numerical analysis, we take tti k = 820 MeV for illustration. 
For the coupling constants c<j and c m , we use [HJ 

\c d \ = 3.2 x 10" 2 GeV, |c m | = 4.2 x 10" 2 GeV, c d c m > 0. (31) 

Although there is no empirical value of c K , we may get it by comparing the ci's obtained in different ways in Eq. (|22[) 
and Eq. @ 

\c K \ = 1.9, c K c m > 0. (32) 
Moreover the coupling constants should obey the nonet relations in the large N c limit, 

C - C C 

Cd — C m = — ^=C m , Cfj — — j=C K: £ = zbl. (33) 

In this way, we get the LECs at 0(e 2 ), 

c = 0.10 GeV" 1 , ci = 0.12 GeV -1 , c 2 = -0.30 GeV" 1 , c 3 = 0.42 GeV" 1 . (34) 

The resonance saturation method may bring large uncertainty in the determination of the LECs at the third order. 
We take the value of n r in Ref [l8[ 

K r = -0.33 GeV" 2 , (35) 

which is obtained by fitting the lattice QCD result s I37J . We simply assume the other tree diagram contributions at 
0(e 3 ) arc small and neglect them as done in Refs. fl2i. |38|. 

V. NUMERICAL RESULTS AND DISCUSSIONS 

We show the numerical results of the T-matrices order by order and the scattering lengths with HM^PT in Table 
Hf. The positive real parts of d^ffi , Q^kd*' a ;?£)*> a i?D*' a rjD^ an( i a ^D* indicate that the interactions are attractive 



4 The numerical values in Ref. flSII also need a few corrections. The corrections are given here in the form of {0(p 2 ), Total, Scattering 
lengths} T . In Table I, they are {-14.2, 1.1,0.04} (0) , {6.4, 7.5 + 5.5i, 0.23 + 0.17i} T£) , and {-6.7, -6.2 + ll.li, -0.19 + 0.35i} Td . 

T DK 

In Table II, they are {-14.2, 0.5, 0.02} T ( 0) , {6.4, 7.3 + 5.5i, 0.26 + 0.20i} Tj5i , and {-6.9, -6.6 + H.lt, -0.24 + 0.35i} Tfl • In Table 

BK 

III, they are {-15.1, -0.6, -0.02} (0) ,{6.1,7.3 + 5.51,0.22 + 0.171}^ , and {-7.4, -7.0 + ll.lt, -0.22 + 0.35i}t d • To get a positive 

T DK 

a D s r] an d nearly vanishing ajy glT , one requires Ci > 0.8 GeV -1 and Co < 4.2 GeV . 
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TABLE I: The threshold T-matrices for the pseudoscalar meson and D* meson scattering order by order in units of fm with 
HMxPT. 





0(0 


0(e 2 ) 


loop 


tree 


total 


Total 


Scattering length 


J 7rD* 


-3.2 


0.5 


-l.-0.0096i 


0.17 


-0.88-0.0096J 


-3.6-0.00962 


-0.13-0.00036i 


T (l/2) 
x -kD* 


0.0 


U.O 


0.53-0.0096i 


-0.33 


0.19-0.009oi 


7.1-0.00961 


0.27-0.00036i 


T (l) 
tD* 





0.09 


-1.1 





-1.1 


-1 


-0.039 


T (o)' 

KD* 


1 ^ 




1 1 n nnm r; 


Q 8 


1 i n nnm 

l.l'U.UUUlUf 


94 n nnm r; 

il. U.UUUlUt 


n o y 1 n~ 6 »' 

U. 1 U J.i A 1U i 


T (l) 





0.75 


-1.5+5.6i 





-1.5+5.6i 


-0.7+5.6; 


-0.022+0. 18i 


T (l/2) 
± KD* 


-7.6 


4.1 


-5.9 


4.9 


-0.98 


-4.5 


-0.14 


T (l) 


-7.6 


4.1 


-7.4-0.000054; 


4.9 


-2.5-0.000054; 


-5.9-0.000054; 


-0.19-1.7 x 10" 6 i 


T fc? * 

JtD! 


7.6 


-2.6 


8.8+0.00016; 


-4.9 


3.9+0.00016i 


8.9+0.00016; 


0.29+5.2 x 10~ 6 i 


7.6 


4.1 


4.+8.3i 


-4.9 


-0.86+8.3; 


ll.+8.3i 


0.35+0.27; 


T (l/2) 

± r 1 D* 





1.2 


0.46+3.; 





0.46+3.i 


1.7+3.i 


0.051+0.094; 


T (o) 





5.8 


0.0036+6.1; 





0.0036+6.H 


5.8+6.1; 


0.18+0.19i 



TABLE II: Comparison of the T-matrices from the loop diagrams for the pseudoscalar meson and D* meson scattering between 
HMxPT and IR. 





Intermediate state: Z) meson 


Intermediate state: D* meson 


Loop: 


total 




HMxPT 




IR 




HMxPT 


IR 


HMxPT 


IR 




5 = 142 MeV 


8 ->■ 


8 = 142 MeV 


5^0 


8 = 142 MeV 


8 = 142 MeV 


8 = 142 MeV 


8 = 142 MeV 


T (3/2) 


-0.053-0.0096; 


0.014 


-0.043-0.0076; 


0.0045 


-0.99 


-0.84 


-l.-0.0096i 


-0.88-0.0076i 


T (l/2) 
A nD* 


-0.053-0.0096; 


0.014 


-0.05-0.0093i 


0.031 


0.58 


0.3 


0.53-0.0096; 


0.25-0.0093; 


T W 
ttD% 


-0.043 


-0.046 


-0.03 


-0.04 


-1.1 


-0.88 


-1.1 


-0.91 


rpW 

1 KD* 


0.69-0.00016; 


0.93 


0.46-0.00015; 


0.76 


10. 


7.5 


11.-0.00016; 


8.-0.00015i 


T w 

KD* 


-0.076+0.000054; 


-0.14 


-0.014+0.0019i 


-0.18 


-1.4+5. 6i 


-3.6+2. 9i 


-1.5+5.6; 


-3.6+2.9; 


T (l/2) 
KD* 


0.46 


0.51 


0.29 


0.43 


-6.3 


-13. 


-5.9 


-13. 


T (l) 
KD* 


0.31-0.000054; 


0.4 


0.17-0.00099; 


0.38 


-7.7 


-14. 


-7.4-0.000054i 


-14.-0.00099; 


T (o) 

KD* 


-0.46+0.00016i 


-0.68 


-0.33-0.0027; 


-0.51 


9.3 


11. 


8.8+0.00016i 


11.-0.0027; 


T (l/2) 
KD* 


0.46 


0.51 


0.31 


0.42 


3.6+8.3; 


1.1+4.4; 


4.+8.3i 


1.5+4.4; 


T (i/2) 

± T}D* 


0.14+0.0021; 


0.16 


0.088+0.0018; 


0.13 


0.32+3.; 


-2.9+1.6i 


0.46+3i 


-2.8+1.6; 


T (0) 
r,D* 


-0.022 


0.013 


-0.015 


0.016 


0.025+6.H 


-6.2+3.2; 


0.0036+6.1; 


-6.2+3.2; 



for these channels. From Table U we see that the chiral expansion of the pion channels converges well. The loop 
diagrams contribute largely to the kaon and eta channels due to the large mass of kaon and eta. But luckily they are 
cancelled by the tree diagram at 0(e 3 ), which makes the whole result convergent. 

We compare the loop contribution between the HMxPT and IR scheme in Table |TTJ For both cases the dominant 
loop contributions are those with the intermediate state D* meson. The numerical results are similar in the pion- 
scattering channels with these two different schemes. But the results differ greatly in the eta scattering channels. 

The mass difference S affects our results only through the intermediate D meson in the loop diagrams. For 
comparison, we list the ^-dependent part of the T-matrices in Table [TT] when 5 — > and S = 142 MeV. We notice 
that the correction from the heavy quark symmetry breaking in the loop diagrams with the D meson intermediate 
state is small in the KD* channels. However, such a correction is significant in the 7rZ?* /=3 / 2 , nD* I=1 / 2 and riD* s I=0 
channels. 

From Table UH the IR scheme lowers the loop contribution in the channels >^i?» ^^d*^kd*' an< ^ ^Rd* ' 

The T-matrices of the 4>D scattering in the nonrelativistic xPT were compared with those in the relativistic X-PT 
with the extended-on- mass-shell renormalization schemes in Ref. 20]. The relativistic effect would also lower the 
loop contribution in some channels such as T^ 2 \t^ 2 \ 

The resonance D s i(2460) couples to the D*K strongly. Its role is similar to that of A(1232) in the case of the pion 
nucleon scattering. Z? s i(2460) is very close to the D*K threshold and may be quite important for the D*K T-matrix. 
In contrast, the non-strange P-wave axial-vector D meson lies well above the D*tt threshold. Its contribution is less 
important. In this work we have tried to include some of the corrections from the P-wave axial- vector D meson 
through the LEC C3. We expect that the results would be improved particularly for the D*K channel if D s i(2460) is 



9 



TABLE III: The threshold T-matrices for the pseudoscalar meson and B* meson scattering order by order in units of fm with 
HMxPT. 







0(e 2 ) 


loop 


0(e d ) 
tree 


total 


Total 


Scattering length 


1 -kB* 


-3.2 


0.44 




1 7 


u.oo 


o.u 


-0.14 




u . o 


0.44 


U.UO 


u.oo 


0.24 


7 1 


u.^o 


7rB* 





0.063 


-1.1 





-1.1 


-1 


-0.04 


T (0) 


15 


6.9 


11 


-9.8 


0.72 


23 


0.83 


KB* 





0.53 


-1.4+5.6i 





-1.4+5.6? 


-0.88+5.64 


-0.032+0.2i 


T (l/2) 

Kb: 

T (l) 

KB* 


-7.6 


3.7 


-6.2 


4.9 


-1.3 


-5.2 


-0.19 


-7.6 


3.7 


-7.6 


4.9 


-2.7 


-6.6 


-0.24 


T (0) 
KB* 


7.6 


-2.6 


9.1 


-4.9 


4.2 


9.2 


0.34 


T (l/2) 


7.6 


3.7 


3.7+8.3i 


-4.9 


-1.2+8.3i 


10.+8.3i 


0.37+0.3i 


T (l/s) 





1 


0.37+3i 





0.37+3i 


1.4+3.i 


0.049+0.11* 


T (o) 

Vs.: 





5.2 


0.022+6. li 





0.022+6.H 


5.2+6. li 


0.19+0.22i 



included as an explicit degree of freedom as A(1232) in heavy baryon chiral perturbation theory [39-4 lj. 

We can also study the pseudoscalar meson and B* meson scattering with P* and P representing (B~ , B° , B®) and 
(B*~ , B*°, B*°), respectively, in Eq. The situation is very similar to the D* case. With the following different 
constants 

c = 0.08 GeV" 1 , a = 0.10 GeV" 1 , c 2 = -0.25 GeV" 1 , c 3 = 0.44 GeV~\ 

K r = -0.33 GeV~ 2 , Af = 5323 MeV, S = 46 MeV, (36) 



and g — 0.52 42j, we can also analyze the interaction of the pseudoscalar meson and B* meson. The numerical results 
with HMxPT are shown in Table HID They are only slightly different from those of the pseudoscalar meson and D* 
meson scattering. 



As mentioned in Ref 18|, it is easy to get pseudoscalar meson and heavy antimeson scattering lengths with the 
C-parity transformation, 

T (I) _rp{I) rp(I) _ T (I) rp{I) _ rp{I) (n 7 \ 

HK HK' HK HK ' 1 H^/r, ~ 1 Htt/t;' ^ ' > 

where / is the total isospin and H(H) denotes the heavy meson(antimeson). 

In short, we have investigated the pseudoscalar meson and D* meson scattering lengths to 0(e 3 ) with HMxPT and 
IR methods. The chiral expansion in the ttD* channels converges well. We hope our present study may be helpful 
to the possible extrapolation in the future lattice simulation of the light meson and heavy meson scattering. Our 
results may also be useful to the phenomenological investigation of the possible molecular states composed of one 
heavy meson and one light meson. 
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Appendix A: The uncertainty of LECs with resonance saturation method 

Other resonances with the same quantum numbers as a, k, and D sl (2460) will also contribute to Cj's. Generally, 
the heavier the resonance is, the less its contribution to the LECs. Here we will check the uncertainty which other 
resonances would cause. 
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From Eq. (1291) . we see that a does not contribute to c\ and C3. For the LECs Co and c 2 , the contribution JC from 
the a and k scales as 

£4 = 4 = 0-3- (AD 

In principle the LECs would absorb the contributions from the /o(1370) singlet and the K$ octet (iiTg (1430), ao(1450), 
/o(1500)) and other heavier resonances similarly. Similar to Eq. (|A1[) . 

|£/o(1370)| = m /o(1370) = Q g ^ A2 ^ 



In the channel /o(1500) — > 7777 and /o(1500) — > KK, the decay momentum is 516 and 568 MeV respectively. We 
can estimate the magnitude of c<i,k* and c m ,K * using the experimental partial width 

r(/ (1500) -> rjq) = 109±£ x (5.1±g;|)% MeV, r(/ (1500) -> = 109±£ x (8.6±J;g)% MeV, (A3) 
and similar effective Lagrangians as in Eq. ([21))) . One gets 



c rf , jf- 1 - 1.9 x 10 2 GcV, |c m ,K* I - 3.0 x 10 2 GcV, c rf> K* c m> K* < 0; or 

Cd,Ks\ ~ L3 X 10_2 GeV ' l Cm '^o* I ~ L5 X 10_2 GeV ' Cd > K a ° m > K a > °" 



(A4) 



Therefore, 



fel _ i^oi ^^L < = M „ 3 (A5) 



I^kI |Cd/ r 



In the last step we have assumed that the coupling ck* in the D*D*Kq vertex is of the same order as c K . Now we 
have ~ 0.3|/C K | and |^-/ (i370) | ~ 0-1 I n other words, the Cjs' correction from the heavier resonances is 

roughly 30%. 

In principle, heavy vector resonances D\ (2420) (D s i (2536)) would also give corrections to C3, 

\Gi + G 2 \ 2 D (2420) M D * 
AC3 = ~2(M 2 — M 2 ) ' (A6) 

By fitting r(.D s i(2536) -> D*K) < T{D sl (2536)) < 2.3 MeV, one finds \G X + G 2 |i» 1 (2420) < 0.16. So |Ac 3 | < 
0.03 GeV" 1 , which is less than 10% of c 3 in Eq. 

Moreover, the variation of the mass of the k octet from 658 MeV to 985 MeV would introduce the 20% uncertainty 
to Co and c 2 and 60% uncertainty to c\ and C3. In short, the determination of q's in Eq. (|34[) are not accurate. But 
their sign and order of magnitude should be reliable. 



Appendix B: Some functions and constants in the infrared scheme 

We perform the tensor decomposition of the IR integrals as follows, 
1 f d d k {l,h> i ,k> i k"} 



i J j (2Tr) d (k 2 -m 2 + ie) [(p - k) 2 - M 2 + ie] 

{l^(p 2 ,m 2 ,M 2 ), p^I^(p 2 ,m 2 ,M 2 ), g^I^{p 2 ,m 2 ,M 2 )+p^p u I^{p 2 ,m 2 ,M 2 )} , (Bl) 
d d k {l,k ti ,k' 1 k' / ,k> i k v kP} 



i 7/ (27r) d (k 2 -m 2 + ie)(k 2 - m 2 . + ie) [(p - k) 2 - M 2 + ie] 
= |F (p 2 ,m 2 ,m 2 ,M 2 ), p^F^p 2 , m 2 , m 2 , M 2 ), g^F 2 {p 2 1 m\, m 2 , M 2 ) + p^F 3 (p 2 , m 2 , to 2 , M 2 ), 

p>V^4(p 2 ,m 2 ,™ 2 , Af 2 ) + (0"V + 5 M V + 9 L/p P^)I r 5(p 2 , ml, m 2 , M 2 )} , (B2) 

where M is the mass of the heavy meson and to, TOi, 7712 are the masses of the light pseudoscalar mesons. The Lorentz 
invariant coefficient can be written as po| . 

I^(p 2 ,m 2 ,M 2 ) = - y - - 2 L~^~ 2 21og--l 

p* 32n^ p z A 
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a 

x 



8tt 2 w 2 



'-VT~ fFcos-^-^) |0|<i 
VH^iog (VC^) 2 - 1 - n<-i, (B3) 

VfF-T {.tt + log (s±s - ^ffi+iZ^-T) } 



n > 1 



where 



* = f ~" - -W ^VTT2^, L^^+lh,-,-^)}. (B4) 
The other coefficients are 

/M^m 2 ^ 2 ) = ^ " ^ 2 2 +m V (p 2 ,m 2 ,M 2 ) + + _i_ log 

j( 2 V,m 2 ,M 2 ) = ^ (m 2 j(°)(p 2 , m 2 , M 2 ) - ^ ~ ^ + I^fr 2 , m 2 , Af 2 )) , 
/( 3 V,m 2 ,M 2 ) = P 2 ~ Af 2 + ™ 2 J( i) ^ ro a M2) _ ]_ J(2) (/> m2> M 2 }> 

— 

^ 4 (p 2 ,TO 2 ,m 2 ,Af 2 ) = -L (^+2/(3)^^^) + ^+2^2 _ M 2 + to 2 )F 3 ( P 2 , m 2 , to 2 , Af 2 ) 

-21 {1 \p 2 , to 2 , Af 2 ) - 2to 2 F 1 ( P 2 , to 2 , m 2 , Af 2 )) , 
F 5 (p 2 ,TO 2 ,m 2 ,Af 2 ) = ^f( 3 > (p 2 , m 2 , Af 2 ) + P ' ~ ^ + ^ F 3 (p 2 ,m 2 ,m 2 , M 2 ) - ^f 4 (p 2 , m 2 , to 2 , Af 2 ). (B5) 
With the definitions, 

. g 2 {2A/ 2 F 5 (M 2 ,,m 2 ,TO 2 ,Af 2 )+4M 2 ,f 5 (M 2 ,,™ 2 ,m 2 ,Af 2 ,)} 
^i(to) = — " ji L > 

h 2 (m) ee j^^f( 2 V,TO 2 ,Af 2 )U^ 



h 3 (m) = ^^^(M^.m 2 ^ 2 )!^^^^^ 2 ^^,™ 2 ^ 2 )!^., 



(B6) 



we can show the functions used for IR in Eq. (|T71) 
£ 



2 

= 7j {4M £ ,.F 2 (M D « , to, to, Af D .) + 2M D F 2 (M D , , TO, TO, Af£))} r , 



q 2 

J = p{4M D ,F2{M D ,,rn v ,m v ,M D ,) + 2M D F 2 {M D ,,rn v ,rn v ,M D )} r , 

V(m, u) = log - - {2A/Z3- w 2 /(°) (w + M D , , to, Af D , ) + 4Af 2 . wl« {uj + M D . , to, M d ,) 

+2M d * I (2) (w + Mr,* , m, M D . ) + 2Af£. f (3) (w + M D * , m, M D . )} , 

Vi = to^ (-^/ii(mRr) - h^m^) - ^h 2 (m^) ~ ^h 2 (m K ) - -^h 2 (m v ) + ^h 3 (m n ) - -^h 3 {m v ) 

( 2 2 12 1 

V 2 = m K <\-hi(m K ) - -h 2 (m K ) - -h 2 (m v ) + -h 3 (m K ) - -h 3 (mr,)j , 

V3 = % ^-hi(m^) - ^h 2 (m n ) - ^h 2 (m K ) - -^h 2 (m r] ) + ^h 3 (m„) - ^h 3 {m K ) - -^h 3 (m n )^ , 

(3 13 11 

V 4 = m K |-/i 2 (m ff ) - h 2 (m K ) - -h 2 (m v ) + -h 3 (m^) - h 3 (m K ) - -h 3 (m v )j . (B7) 
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The {A"} r in Eq. (|B7I) represents lim X after removing the terms proportional to L. 

d—>4 
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